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We derive spin-orbit coupling effects on the gravitational field and equations of motion of compact 
binaries in the 2.5 post-Newtonian approximation to general relativity, one PN order beyond where 
spin effects first appear. Our method is based on that of Blanchet, Faye, and Ponsot, who use a post- 
Newtonian metric valid for general (continuous) fluids and represent pointlike compact objects with 
a 5-function stress-energy tensor, regularizing divergent terms by taking the Hadamard finite part. 
To obtain post-Newtonian spin effects, we use a different 5-function stress-energy tensor introduced 
by Bailey and Israel. In a future paper we will use the 2.5PN equations of motion for spinning 
bodies to derive the gravitational-wave luminosity and phase evolution of binary inspirals, which 
will be useful in constructing matched filters for signal analysis. The gravitational field derived here 
may help in posing initial data for numerical evolutions of binary black hole mergers. 
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I. INTRODUCTION 



Gravitational waves from coalescing compact binaries are the most promising candidate for detection by near- 
future, ground-based laser interferometers such as LIGO, VIRGO, GEO600, and TAMA [jjj. Detection of signals 
and estimation of signal parameters in noisy data require detailed modeling of the inspiral phase of the waveform 
H and — depending on the mass of the system — some knowledge of the merger phase ||. The matched filtering 
techniques used to analyze the inspiral waveform require that the model, or template, match the real signal to less 
than one cycle out of the roughly 10 3 -10 4 in the detectable band. This precision requires high-order post-Newtonian 
calculations of the equations of motion and gravitational-wave luminosity. The late inspiral and merger stages of 
coalescence are not amenable to post-Newtonian approximation methods and require numerical evolution of binary 
black-hole spacetimes. However, a post-Newtonian approximation of the gravitational field may help in posing initial 
data for such numerical evolutions B. 

The standard method for modeling inspiral waveforms from binaries of arbitrary mass ratio is the post-Newtonian 
expansion in powers of the binary's orbital velocity v/c and gravitational potential GM/rc 2 ~ (v/c) 2 . There are 
two nearly equivalent approaches to calculating gravitational waves at high Post-Newtonian order: one developed 
by Blanchet, Damour, and Iyer (BDI) |3]-[7| and one by Will and Wiseman || based on previous work by Epstein, 
Wagoner, and Will ||. The gravitational waveforms and luminosity are expanded in time derivatives of radiative 
multipoles, which are expressed as integrals of the matter source and gravitational field. The radiative multipoles 
are combined with the equations of motion to obtain explicit expressions in terms of the source masses, positions, 
and velocities, which can be converted to gauge- invariant frequencies observed at infinity. Under the assumption that 
the bodies can be treated as pointlike particles characterized only by their masses, the matter source (stress-energy 
tensor) is given as a (5-function. In combination with a regularization scheme for infinite self-field effects, this source 
greatly simplifies the field integrals compared to a detailed fluid body calculation. 

It is standard to count post-Newtonian (PN) orders in powers of (v/c) 2 beyond the Newtonian result for the equation 
of motion and the quadrupole formula for the gravitational waves. For the case of nonspinning bodies, the equations 
of motion have been known to 2.5PN order for decades |Tc|-|l3||, and the gravitational- wave luminosity also has been 
evaluated recently to 2.5PN order m. At the moment the 3PN calculations of the equations of motion for nonspinning 
bodies are nearing completion (see |14] , |i~5t and references therein). 

The post-Newtonian expansion of gravitational waves from spinning bodies has not been carried as far as for 
nonspinning bodies — almost exclusively the literature contains calculations of leading-order effects. There are several 
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kinds of spin effects: Spin causes precession of the orbital plane of a binary, changes the orbital frequency, affects 
the gravitational-wave luminosity, and modifies the amplitudes of the gravitational waveforms. Spin effects can be 
further divided into spin-orbit coupling (involving a single spin) and spin-spin interactions. f] The leading-order terms 
in the equations of motion were derived in the 1950s by Papapetrou flq] , and in precession in the 1970s by Barker and 
O'Connell These results have been re-derived by several methods [ ]lT| , |l8| -p2"[ . Kidder, Will, and Wiseman |^3| in 
the 1990s evaluated the leading-order (1.5PN spin-orbit and 2PN spin-spin) terms in the gravitational-wave luminosity 
(not including precession). Kidder [ p4[ completed that work with evaluation of the lea ding -order (1PN spin-orbit and 
1.5PN spin-spin) terms in the waveform amplitudes. He, along with Apostolatos et al. p5fl , also considered the effects 
of precession on the waveforms. All of these calculations were carried out only to leading post-Newtonian order in each 
spin effect, due in large part to the unwieldiness of calculations representing spinning bodies with a fluid stress-energy 
tensor of nonzero compact support. 



However, it is possible to make a "spinning particle" approximation to simplify the calculations. Mino et al. [|26 27 
used a 5-function stress-energy tensor, a compact reformulation of the work of Dixon [ p8| , to represent a spinning test 
mass in orbit around a massive spinning black hole to high post-Newtonian order in the limit of small mass ratio. Cho 
p9| used a similar stress-energy tensor due to Bailey and Israel J3(J to re-derive the lowest-order results of Kidder, 
Will, and Wiseman |23|. Recently we j3l[] made the first calculation in full (arbitrary mass ratio) post- Newtonian 
theory of a non-leading order effect, a 2PN spin-orbit contribution to the waveform amplitudes, using the stress-energy 
tensor of Mino et al. pq | . Our goal is to eventually obtain the lowest-order unknown spin effect in the phase evolution 
of the waveforms, a 2.5PN spin-orbit term. The new term requires evaluation of the radiative multipoles at relative 
1PN order, i.e. 1PN beyond the leading-order 1.5PN effect, a task which can be done in a straightforward manner 
using the techniques of Ref. |3l| ] . However there is also a contribution from the 2.5PN spin-orbit terms in the equations 
of motion, which we must first calculate. 

Recently Blanchet, Faye and Ponsot (BFP) |3^] derived the equations of motion for nonspinning objects using 
a new approach, based on the post-Newtonian fluid metric used in wave generation calculations combined with a 
(^-function source and the use of the Hadamard finite part |l3|,|33| to regularize the resulting divergent integrals. In 
this paper we generalize the BFP calculation to include spin, thereby obtaining the missing 2.5PN spin-orbit term 
in the equat ions of motion. This term, together with the well-known 2PN spin-spin term and the 2PN quadrupole 



term []34 35 1 completes the equations of motion of two finite bodies to 2.5PN order. The calculation is similar to our 
previous work pl| , except that we now use the stress-energy tensor of Bailey and Israel p0| due to its advantages in 
deriving the regularized equations of motion and precession. 

We have made some effort to check that our spinning particle approximation is actually a good model of a spinning 
compact object. Of course, any <5-function approximation of an extended body yields nonsensical results near the 
body, but we only require a good approximation to the field some distance outside the body and to quantities obtained 
by volume integrals of the field or its source. (The use of our field results as initial data for numerical binary black 
hole simulations must be supplemented by an additional prescription near the bodies, such as matched asymptotic 
expansion jij.) Following the lead of Ref. p3[] , we give some consistency arguments. First, we show that the one-body 
limit of our result for the gravitational field reproduces the Kerr metric up to the post-Newtonian order considered. 
Second, we verify that the equations of motion and precession we derive from a regularized version of the stress-energy 
conservation law satisfy the harmonic gauge condition assumed at the beginning. Third, in our approach we reproduce 
the well-known 1.5PN spin-orbit terms in the equations of motion. Finally, we verify that the new 2.5PN terms we 
derive are Lorentz invariant and reduce (in the test mass limit) to the results of black-hole perturbation theory. 

This paper is organized as follows. In Sec. II, we describe our spinning particle approximation, including the stress- 
energy tensor, regularization scheme, and derivation of the regularized equations of motion and precession. In Sec. 
Ill, we use the post-Newtonian expansion to write needed quantities in terms of a set of post-Newtonian potentials. 
Section IV is a description of the calculation of the post-Newtonian potentials in terms of properties of the bodies. 
The (lengthy) results for the potentials are relegated to Appendix A. In Sec. V, we present the new 2.5PN terms in the 
equations of motion. In Sec. VI, we summarize our results and describe problems of future interest. In Appendix B, 
we check our result for the circular-orbit frequency by comparing it to results for spinning test particle in orbit around 
a Kerr black hole. In Appendix C we provide a brief summary of the post-Newtonian check for Lorentz invariance. 

Throughout this paper, we use units such that Newton's gravitational constant and the speed of light equal unity. 
The brackets () and [] on tensor indices indicate symmetrization and antisymmetrization respectively: 



At higher post-Newtonian orders it is also possible to have three-body, three-spin interactions; but these have not been 
considered yet. 
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Greek indices run from to 3, and Latin indices from 1 to 3. 

II. SPINNING PARTICLE APPROXIMATION 

In this section we describe how we simplify fluid-body calculations by approximating spinning compact objects as 
pointlike particles. We use the Bailey-Israel (^-function stress-energy tensor and the Hadamard finite-part regulariza- 
tion scheme for the gravitational field of a point source. Conservation of the stress-energy tensor and vanishing of 
its antisymmetric part yield the Papapetrou equations of motion and precession in a form where the regularization is 
unambiguous. 

A. Stress-energy tensor 

The particle-like stress-energy tensor of Bailey and Israel |3C| l can be divided into monopole, spin, and spin anti- 
symmetric parts: 

rpaff _ rpUfi , rpUfi , rriQ/3 Icy 1 \ 

1 ~~ (M) (S) ' (SA) ' V- k ) 

These parts are expressed as integrals over the trajectory y(r) of the particle as 

= / dr^(r)t^(r) f(4)( ^- y(T)) (2.2a) 



(2.2b) 



T" s %{x) = - iv 7 / dr^rMr) ^^ - (2.2c) 

Here w M = dy^/dr, p^{t) is the linear momentum of the particle, and S a ^(T) is an antisymmetric tensor representing 
spin angular momentum. However, the trajectory parameter r is generally not quite the proper time of the body (see 
below). To re pre sent multiple bodies, we label each particle and its associated quantities by A (y 1 ^, etc.) and sum all 



parts of Eq. (2.2) over A. 

When treating spinning compact bodies as point particles, we introduce spurious degrees of freedom which must 
be fixed. A finite spinning body has a center of mass (which must be defined carefully in post-Newtonian theory, see 
e.g. Ref. p3[). However, in the limit as the body is shrunk to a point particle, the trajectory of the particle does not 
necessarily correspond to the worldline of the finite spinning body's center of mass. The ambiguity in precisely what 
the trajectory represents (or cquivalcntly in defining the center of mass) is fixed by imposing a spin supplementary 
condition. We choose 

S^pu = 0, (2.3) 

but note that other choices are common in the literature (see the discussion by Kidder p4||). These other choices of 
spin supplementary condition lead to equations of motion that appear different but are physically identical, a common 
source of confusion. The ambiguity in the definition of the spin also implies that r is not precisely the proper time 
along the body worldline (although r reduces to the proper time in the Newtonian limit). Thus u M (r), the tangent 
vector to the trajectory y^(r), is not precisely the four- velocity of the body. We define 

u^^p^/m, (2.4) 

where m is the Schwarzschild mass of the body. Holding the linear momentum constant along the worldline implies 

g^u v = -1. (2.5) 

We also introduce a spin vector four- vector S 1 *. The spin supplementary condition (|2.3j) allows the natural definition 
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s 



a/3 _ 



-<] 



(2.6) 



where e^ vp(y is the totally antisymmetric symbol with g 0123 = 1. Even after demanding that S^S^ be conserved along 
the trajectory and that S l reduce to the Newtonian spin in the appropriate limit, S° is still undefined. We fix this 
remaining degree of freedom by imposing the condition 



S^p^ = 0. 

Note that the literature also contains different conventions for this condition (see e.g. |20|]). 



(2.7) 



B. Regularization procedure 

In deriving the equations of motion, we need to evaluate the gravitational field and associated quantities at the 
locations of the bodies. However, as usual when dealing with pointlike sources in field theory, the quantities we need 
diverge if we naively integrate the (unphysical) (^-function stress-energy tensor. Thus we must augment our method 
with a regularization procedure to remove the infinite self-field effects while preserving the correct physical terms. In 
this paper, we adopt the Hadamard finite part p^ , p2| , p3| , p6[ . 

In our calculations we encounter a class of functions F which admit, when the field point x approaches one of the 
source points {r\ = \x — y\ | — ► 0), an expansion of the form 



F(x)= ^/fc("i)+0(ri), 



(2.8) 



-fe o <fe<0 



where k is a set of integers and ri\ = (x — y\)jr\. We define the effective value of the function F at the position of 
body 1 as the Hadamard finite part, 2 

(J)i=%,»i) = / ^^/o(m,yi). (2.9) 

This effective value also holds for integrals of F with a three-dimensional (5-function, 

J d 3 xF(x, y x )5{x - Vl ) = (F), . (2.10) 

We introduce a symbol Qa to denote the Hadamard finite part at x — tja of the function within parentheses (). 
When we calculate the post-Newtonian potentials, encounter integrals of the form 

J d 3 xF(x)T a(3 . (2. II) 

Each part of this integral can be evaluated by means of Hadamard's regularization as 

d 3 xF(x)T$(x) 



(2.12a) 



E 



s A y -p= 



dt sTv 



oK,/ 3 ) 



-9 J A 



1 ~f5 A V A 



y(a r /3) x\ F 



(2.12b) 



2 A more sophisticated treatment of the Hadamard finite part is used in some higher-order calculations [^i). However, this 
version, which is based on that of BFP, is sufficient for our purposes. 
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d*xF(x)T«%(x) 



2 V 



Si"' 



'A V A 



F 



9 J A 



A V A 



F 



(2.12c) 



C. Equations of motion and precession 



The equations of motion and precession can be derived from the regularized stress-energy tensor. However, the 
regularization does not commute with some derivatives and multiplications, and moreover our version of it takes 
place on slices of constant t rather than in the rest frame of each body. Therefore at each step we must make some 
consistency checks. 

We consider the following form of the conservation of the stress-energy tensor, 



-gT 



a/3\ 



+ P 



-gT^ = 0. 



(2.13) 



By using the Bailey-Israel stress-energy tensor (2.2) and integrating over a three- volume Da containing only body A, 
we obtain the equations of motion as 



1 



dt 



Pa + {^) a P>a + ^ (R a ^)y A ST = 0, 



(2.14) 



where v A = u A /u A . (If we used Dixon's stress-energy tensor as before [[51J, we would have both the equations of 
motion and of precession entangled in this expression, and it would be difficult to separate them.) By performing the 
same integral with the free index covariant, we obtain 



d_ 

dt 



PA a - {^ av ) A pA^v A + - (R a p^) A v A S^ = 0. 



(2.15) 



It can be shown by direct calculation that (2.14) and ( 2.15| ) are equivalent up to 1PN order for non-spinning terms 
and 2.5PN order for spinni ng te rms, which is a useful consistency check of our calculations. Direct calculations also 
show that Eqs. (2.14) and (2.15) are Lorentz invariant up to 1PN order for non-spinning terms and 2.5PN order for 
spin-orbit terms. 

To derive the regularized equations of precession from the Bailey-Israel stress-energy tensor, we demand that 



d 6 x^=g~T [al3] = 



(2.16) 



[Strictly speaking, we require that h\ a ^ vanish, but to the order we need this condition reduces to Eq. ( 2. 16 ).] It is 
straightforward to derive the spin precession equations from this volume integral as 



f t sf + (rvh Sfv A + (T^JTva + a?M = o. 

In the same way, we can also derive the spin precession equation with covariant indices as 



d_ 

di 



(s AaP ) A + {T^) A SA^vl+(l 



Pi 



SAa^V A 



+1pA[aVA0] = 0. 



As we shall see in later sections, ( 2.1 7| ) and ( 2.1§| ) are equivalent, at least to the order we need. 
Using Eqs. ( |2.14p and fl2.17| ), we can also derive the relation 



Pa ~ mAU A = 



2mA 



{Rfjupa)^ S A u A S p A . 



(2.17) 



(2.18) 



(2.19) 
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Equation J2.19 ) tells us that if we consider only spin-orbit effects we can neglect the difference between and 
Therefore in the rest of this paper we set 



and set r to be the proper time of the body. 



dT 



(2.20) 



III. THE POST-NEWTONIAN EXPANSION 

We introduce a parameter e of the order of a characteristic velocity v. Since we consider bound systems, we can use 
the virial theorem to put v 2 — 0{M/r) for the total mass M and (harmonic coordinate) separation r. To isolate spin 
effects, we also introduce a parameter x which is also dimensionless and (at most) of order unity for compact objects. 
This x is °f order (spin) /(mass) 2 for some combination of the bodies. Terms of 0(x) correspond to spin-orbit effects, 
spin-spin effects are 0(x 2 ), etc. We use the shorthand 0(m,n) to denote terms of order e m and xe™, and 0(m) to 
denote terms simply of order e m or x £m (depending on the context). In this paper we neglect terms of 0(x 2 ) and 
higher. 

We will expand numerous quantities in s and Xi but first let us determine what orders are needed. In Eq. (2.14), 
the Newtonian force is 0(e 4 ); thus the 2.5PN spin-orbit terms are 0(xe 9 )- By expanding the entire expression in 
terms of the metric tensor, we find that we need <?oo to 0(4, 7); gio to 0(5,6); and gij to 0(4, 5). By integrating 

we find 

1 I J ( 



different terms in the stress-energy tensor ( |2.2| ) over 3- volume and substituting S ~ x m " 

0(0,3), 
0(1,2), 



(AO 



1 / 1 (M) 



T«/T^ =0(2,3). 



We can insert this order counting into the metric below. 

The Einstein equations can be written in harmonic coordinates as 



where 



and the harmonic gauge condition is 



EE 



-ggT - vT 



(3.1) 



(3.2) 



(3.3) 



h^„ = 0. 



(3.4) 



Here □ = rj^d^d^, rf v — diag(— 1, 1, 1, 1), and k^ v (h) represents the non-linear terms in the Einstein equations. 
It is convenient to define the densities 
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^00 _|_ rp 




EE T i0 , 




EE T lj , 



(3.5a) 
(3.5b) 
(3.5c) 



and the retarded potentials 



By counting orders with the post-Newtonian metric used in Sec. [II, we see that this difference does not matter until the 



3PN spin-spin terms in the equations of motion and the 2PN spin-spin terms in the equations of precession. 



G 



□ ^ {-4 7 rG((T ij - 5 t3 a kk ) - d t Vd 3 V} , 
a' 1 {-^G(Vo l - 7cx) - 2d k VdiV k 
3. 

i" = n^ 1 {-AkGVctu + 2Vid t diV + Vd?V 



V 

Vi 

Ri 



:d t VdiV 



where 



□^{-47r/(x,t)} = J ^-f( z ,t-\x-z\) 



(3.6a) 
(3.6b) 
(3.6c) 

(3.6d) 
(3.6e) 
(3.7) 



Inserting the order counting of Eq. ( |3.lD , we find that the potentials have the following post-Newtonian orders: 

7 = 0(2,5), 7 = 0(3,4), Wy = 0(4,5), 

A,=0(5,6), 1 = 0(6,7). (3.8) 
To the order we require, the solution to the Einstein equations ( |3.2j ) is given by 

goo = -1 + 2V -2V 2 + 8X + 0(6,8), (3.9a) 

5*o = -47 - 8R, + 0(7, 8), (3.9b) 

9ij = % (1 + 27 + 27 2 ) + iWij + 0(5, 6), (3.9c) 

7^9= l + 27 + 47 2 + 2F7 fefe +0(6,7). (3.9d) 
The harmonic coordinate condition ( |3.4|) reduces to the following identities between the potentials: 



idtVi 



4d( 



\k W kk ) =0(5,6), 



d t (V + \W kk + 27 2 ) + di 
= 0(6,7) 



7 + 2R t + 277 



(3.10a) 



(3.10b) 



(Actually, these are redundant and we can use only the former.) 

The (regularized) post-Newtonian metric (3.9) allows us to relate the relativistic source parameters to the Newtonian 
ones. From here onward, indices on the right-hand sides of expressions are raised and lowered with the Cartesian 
metric Sij and are written up or down merely for convenience. From Eq. (2.5) we obtain 



u° A = l+ \v\ + (V) A + \v\ + b -v\(V) A - 4*4 (7) y 



+ -(7 2 ) A + 0(6,7). 



(3.11) 



From the definition (2.7) of S° we obtain 

S°a = S A [v\ - 4 (Vi) A + iv A (V) A + 0(5)] . (3.12) 

We introduce for conven ienc e the s hor thand S A = e l ^ k S A (where e l ^ k is the three-dimensional antisymmetric symbol), 
and by combining Eqs. (2.3) and (2.6) we obtain 



diQ _ nij 
°A — °A 



A 

"(5)}, 



1 + ^1 



-3(7), 



mi 



a A — ^A 



e ljk v k A (v A SA) 



+m A 0( X e 4 ). 

Here we have also introduced the shorthand for the (Cartesian) scalar product of two vectors (a b) = a % b l 



(3.13a) 



(3.13b) 
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IV. CALCULATION OF THE POTENTIALS 



We now describe the calculation of the potentials to the needed orders. The results are lengthy and are displayed in 
Appendix A. Two of the potentials (V and Vi) are needed to relative 1PN order, i.e. to 0(2) beyond the leading spin 
and non-spin terms, while the rest are only needed to the lowest order at which they appear. To obtain the explicit 
results in terms of the bodies' coordinates, velocities, masses, and spins, we must substitute lower-order results for 
the equations of motion and precession as well as for V and Vi. 

The lower-order quantities we need are easily obtained. The lowest-order potentials are simple to evaluate as 



V = 



V 



mi 
r\ 
mi 



2S^n\v{ 



27T 1 



1^2 + 0(4,7), 
1 ^2 + 0(5,6), 



(4.1a) 
(4.1b) 



where ta = \x — ua\ and n A 

Alternatively they can be calculated by inserting (4.1) and the results of Sec. Ill into Eq. (2.14) or (2.15) as 



(x l — y l A )/rA- The lower-order terms in the equations of mo tion arc w ell known 



dv\ 
~dt 



m 2 
-4(viv 2 ) 



5mi 



4to9 



-(^12^2) 



ri2 

2 



+ ( Vl ) 2 + 2(v2f 




- v{ 2 [4(ni2Vi) 

12 

ri2 

12 (3n k 12 n l 12 -5 hl ) 



V mi mo / ri2 V 



-0(4,5)), 



(4.2) 



with ni2 = \yi — y 2 |/ri 2 and v\2 = vi — v 2 . The le ading order terms in the equations of precession can be calculated 
by the same procedure using Eqs. ( 2.17 ) and ( 2.18 ) as 



^± = H!a {S\{n l2 v l2 ) + 2n\ 2 [(S lV2 ) - (S lVl )} 
at rf 2 

+(Sxn 12 )(vi-2vi) + SxxO{4)}. 



(4.3) 



The equations of motion and precession for body 2 are given by exchanging labels 1 and 2 in Eqs. ( |4.2| ) and ( |l.3] ). 
Note that the effect of spin prece ssion first appears at 1PN order, in the sense that dS/dt — (S v)/r x 0(2). 
The spin precession equations (4.3) are different from the usual form of the geodesic precession pO|,p7| 



dS\ _m 2 t 
— _— £ijkbl e 



:vq (2v\ 



P\ Q 



(4.4) 



This is due to differing definitions of spin. Equation (4.4) holds for spins defined in the local asymptotic rest frame 
of each body such that S\S\ is constant along the trajectory. However, for our definition of the spin vector this does 
not hold true (beyond leading order). The spin definitions are related by 



S{ 



1 



r\2 ' 



Si x 0(4). 



(4.5) 



The compact parts of the potentials (integrals with compact support) are straightforward to evaluate by a retar- 
dation expansion of the integrands of Eqs. ( |3.6a )-(3.6e) and substitution of the lower-order equations of motion (4.2) 
and precession (4.3). We integrate only the symmetric part of the stress-energy tensor. The equations of precession 
(4.3) ensure that T^g^ does not contribute the metric. It can be verified term by term that T"g A ^ does not contribute 
to the individual potentials either, to the order considered here. 

We also encounter what are generally called quadratic terms (integrals of products of two potentials or their 
derivatives) which do not have compact support. We use the distributional derivative 
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2 / 1 \ 3n\n\ - S» Air 

9i i{n) = 71 y*^(*-»0. 

which is sufficient to our desired order fl32| , change time derivatives to spatial derivatives by 



(4.6) 



(4.7) 



where du = d/dy\, use di(l/r{) = —du(l/ri), and move the operator du outside the integral sign. We also make 
frequent use of the identity 



-if 1 



\rir 2 



ln(ri + r 2 + r 12 ) 



(4.8) 



and its derivatives, where A here represents the Laplacian. 



We check the potentials in two ways. First, we verify that they satisfy the identities (3.10b) equivalent to the 



harmonic gauge condition. Next, we verify that the metric reduces to Kerr in the appropriate limit. By using the 
potentials in Appendix A and setting m 2 = 0, v\ = 0, S 2 — 0, and v\ = 0, we obtain the metric 



2mi 2m? ^, 
9oo = -i + — - -nr + °( 6 > 



9oi 



■\n\Sf + ^i n f5 1 w + 0(6,7), 



2mi mf 



9ij = Sij [ 1 + + -± + -if n l X + 0(5, 1 



(4.9a) 
(4.9b) 
(4.9c) 



This is equivalent to the Kerr metric in harmonic coordinates truncated at the indicated post-Newtonian orders. Thus 
we have some evidence that our spinning particle approximation produces physically sensible results. 



V. EQUATIONS OF MOTION AT 2.5PN ORDER 



Using the post-Newtonian potentials derived in the previous section, we derive the spin-orbit interaction term s in 
t he eq uations of motion at 2.5PN order. Direct calculation shows that it does not matter whether we use ( 2.14 ) or 



A. Body-centered form 



We can write the equations of motion in the Newtonian-like form 



V 1 — T % 



The post-Newtonian "linear momentum" is given by 



V\ = m lV \ 
1 



+ l(V 2 ) 1 -4v{(V J ) 1 + 0(6,7) 
The "force" has two components, the "frame force" 

= mi«; (v A (1 +vl - 4V) + 4Vi(l - 2V) + 4Vj (v + 2v{v^ + 8R, + AX, t + HVjVj. 



(5.1) 



(5.2) 



AViWi: 



+v{V d + 4v{ 2V m (1 - 2V) + AR^j] - W k 



2v 



3 A (W jk ,i - 2W» i ,fe)} i + 0(7, 10) 



(5.3) 







and the spin-curvature coupling force 



■pi _ 1711 (pi \ „Pq» v 



mi {S{° (V + VjVj + 4v,j) + S{ j [Vj + W k V [m + VVj + v{V k V k + (V jk - VjV k )\ + S{ 
x Vijil - W) + 4V- (iJ) - ZViVj - v{Vi + 4wf V[ itfc]j ] + S{ k [2V j>ik (l - 2V) + 2W ij}k + AR^ lk 
-8V (i V j)>k + v\ (V t3 - V t V 3 ) + 2v[ (w ijM + W khi ^~\ } + 0(7, 10). 



(5.4) 



Note that our conventions are slightly different from those of BFP. 

By inserting the explicit formulas of the post-Newtonian potentials, given in Appendix A and by performing the 
Hadamard regularization, we obtain the equations of motion in terms of body masses, spins, positions, and velocities. 
The equations of motion of body 1 can be expressed as 

dv{ m 2 r , . \ it, \ 

-7T = — |A)(yi - V2)+Ai{vi - V2,vi,v 2 ) 

at r 12 

+ B iAyi - y^, V! - v 2 , Si, s 2 ) 

+Bi(y l -y 2 ,S 1 ,S 2 ) 
+Bi 5 ( yi - y 2 , v u v 2 , Si, S 2 ) + 0(4, 6)] . 

Here A l and A\ are respectively the Newtonian and 1PN forces independent of spin, given by 

4* — — n i 

A a — "12: 
A\ = n l2 



mi m 2 o o 

1 - v\ - 2v\ + 4( Vl v 2 ) 



5— +4 

ri2 r 12 



:(ni2V 2 y 



+ v[ 2 [4(m 2 ui) - 3(ni 2 u 2 )] 



(5.5) 
(5.6) 

(5.7) 



The leading-order (1.5PN) spin-orbit interaction force is given by 



mi 


+ 2 


S 2 k \ i 
V 12 

m 2 J 


'( 




+ *JL) 




mi 


m 2 1 



L 12) 



(5.8) 



The leading-order (2PN) spin-spin interaction force is given by 



B 2 = — — 

T 12 



~5«i 



(5i5 2 ) (n 12 S 2 ) S 2 (n 12 Sx) 



ni2- 

mim 2 mi m 2 
{n^Sx) (n 12 S 2 ) 



m 2 mi 



mi 



m 2 



(5.9) 



All of these terms are well known and can be found in references such as Damour fl2l| . 

The new term, the 2.5PN (1PN beyond leading order) spin-orbit coupling term, is given by 



B 2 . 5 = 



mir 



{ 14mi 9m 2 \ 



V r 



15 
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23mi 12rn 2 



26mi 18m 2 



15(nv 2 ) 2 - 6(vi2V 2 ) 



2r 
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-\^v\ 2 
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49mi 20to 2 
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mi 



mi 



(5.10) 



B. Center-of-mass form 



Here we express the equations of motion in terms of the bodies' relative coordinate, in the center-of-mass frame. 
We define the mass parameters M = ni\ + m 2 , n = niim 2 /M 2 , and A = (mi — m 2 )/M; and the dimensionless, 
symmetrized spin parameters B 



It is also convenient to introduce the shorthand 



1 (Sj_ 

2 \mj 

1 /Si 



S_ 

m 2 
S 2 



2 V TO i TO I 



2 ; 







+ 4 








\ 








1 mf 





(5.11a) 
(5.11b) 

(5.12a) 
(5.12b) 



The relation between the body coordinates yi, y 2 , and the relative coordinate is 



Vi = x 
V2 = x 



m 2 , 1 



mi 1 



71/ 



71/ 



77A 



71/ 



Mnv x (xa + Axs), 
- Mnv x (xo + Ax a ). 



(5.13) 



The last terms in the above equation are 1.5PN order and contribute to the 2.5PN relative acceleration through 
corrections to A\. From now on we drop the subscript 12 on r, n, and v. 
We define the relative acceleration as 



_ dvi dv 2 
a ~~dt dT' 

and we express the equations of motion as 

a = a N + a PN + a S o + a S s + a PNS o + 0(6, 
The Newtonian acceleration ctjv is given by 

M 

a N = ^n. 



(5.14) 



(5.15) 



(5.16) 



The 1PN acceleration apjy is given by 



M 



apN = -— \ n 



(1 + 3n)v 2 - 2(2 + n)— - %r){nv) 
r 2 



-2(2 - rj) (n v)v} . 
The 1.5PN spin-orbit acceleration a so is given by 



(5.17) 
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M 2 



aso 



{6n(n x v) ■ (xs + Ax<0 



-2v x [(2 - r,) X s + 2A Xa ] 
+6(nv)n x [(1 - i^Xs + &Xa}} 



(5.18) 



The 2PN spin-spin acceleration ass is given by 



M 3 

ass = ~ — 3t? {n [| Xs | 2 - |Xa| 2 - H^Xs) 2 + 5(n X a) 2 ] 
+2[Xs(n 12 Xs) - Xa(nXa)}} ■ 
The new 2.5PN spin-orbit acceleration apjyso (i-e. the post-Newtonian correction to ago) is given by 



(5.19) 



M 2 

a PNSO = ^3 ( XV <j « J ") 



(32 - 15?y + S77 2 ) — + 1577(1 + 2r))(n v) 2 - 677(3 + 2rj)v 2 



M 

(32 + I377) h 15t7(tztj) 2 - I877V 

r 



3(-8 + r?)- 
r 



-3?7(5 + 277) (nv) 2 + Ur/v 2 ] } + Axi 1 \ n 1 (n v) 
-lbT]{nv) 2 + Ur/v 2 ] } + n i n i v k jxf 



M 

-3(8 + 377) — 
r 



(-40 - 1377 - I677 2 ) — - 3077(n v) 2 + 2Arjv 2 



jk 



(-40 - 2I77)— - 3077(n7j) 2 + 24777J 2 
r 

M 

-8*7 [{nxs) + A(nxa)] 677(71 u) (1 + 2?7)(i; x«) - A(i>Xa)] 

r 



6uW(nu) [(2 - 377 - 2t? 2 ) xf + (2 - 7?)Axi fc ] + (n x «)* 

(5.20) 



C. Quasicircular orbits 



When the spins are aligned or anti-aligned with the orbital angular momentum, the spin vectors and the orbital 
angular momentum vector do not precess. In this case, a quasicircular orbit exists. We define Xa = ±|Xa| and 
Xs = ±|Xs|- The signs of Xa & n d Xs are positive when the vectors are aligned with the angular momentum axis 
(to x v) and negative when antialigned. For quasicircular orbits, the equations of motion take the form 



-LU 2 X, 



(5.21) 



where u> is the orbital angular frequency. Using the identity v = tor, we can eliminate v from the equations of motion 
and express to as an expansion in M/r, 

= — { 1 - (3 - 77)— - 2 [(1 + 77) Xs + A Xa ] — 



+^ 2 tfc) 2 - (Xa) 2 ] + [(9 - 377 + 77 2 )X, 



-(9 - 677)A Xa ] 



■0(4,7) 



(5.22) 



In Appendix B, we show that Eq. (5.22) agrees in the limit 77 — + with the motion of a test particle around a Kerr 
black hole to 2.5PN order. 

The new terms in Eq. Q5.22 ) are proportional to M/r 3 (M/r) 5 ^ 2 : The 0(j] 2 ) term and parts of the 0(77) terms have 
not been derived before. Their magnitude indicates that 2.5PN spin-orbit contributions change the circular orbit 
frequency by an amount of less than about two percent near the innermost stable orbit for a system of maximally 
spinning black holes. However, the importance of this effect on the phase evolution of the orbit, and thus on the 
effectiveness of matched filtering for the gravitational-wave signal, cannot be evaluated until the gravitational-wave 
luminosity is known. 
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VI. SUMMARY 



We have calculated the 2.5PN spin-orbit effects on the gravitational field and equations of motion of compact binaries 
by approximating the compact bodies as spinning point particles. Although the formal basis of this approximation 
remains uncertain, we have demonstrated that it is reasonable in that it reproduces some physical behavior of the Kerr 
spacetime to this post-Newtonian order, as well as lower-order, previously known terms in the equations of motion. 

Using the techniques of this paper, it is a straightforward though lengthy task to calculate the equations of precession 
to 2PN order and the gravitational-wave luminosity (needed to construct templates) to 2.5PN order. We plan to do 
this in the future for circular orbits. The 2.5PN gravitational field given by combining our potentials with the spinless 
potentials given by BFP could serve toward providing initial data for the numerical evolution of spinning binary 
black-hole mergers, if augmented by some scheme for approximating the field near and inside the black holes 
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APPENDIX A: EXPLICIT POST-NEWTONIAN POTENTIALS 



In this appendix, we present explicit expressions (in terms only of the masses, spins, positions, and velocities) for 
the post-Newtonian potentials necessary to obtain the equations of motion and the metric away from the bodies. 
We denote spin and monopole parts of quantities with subscr ipt (S) and (M), res pect ively: e.g., Vrs) and Vim)- All 
expressions are fully reduced using the equations of motion (4.2) and precession (4.3). The monopole terms in the 
potentials are given to higher order by BFP, but for convenience we reproduce the pieces needed in our calculation. 

The relevant potentials are 
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where s = n +r 2 +ri 2 - These potentials, inserted into the post-Newtonian metric (3.9), give the spin-orbit interaction 
terms in the gravitational field to 2.5PN order away from the bodies. 

Using the Hadamard finite part, we find the values of the spin parts of the potentials at body 1 to be 
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APPENDIX B: COMPARISON WITH KERR ORBITAL FREQUENCY 



We perform one check of our equations of motion by comparing the orbital frequency (5.22) in the test-mass limit 
to the orbital frequency of a spinning test particle in a circular orbit around a Kerr black hole derived by Tanaka et 
"/• 

The transformation from Boyer-Lindquist coordinates (tsLifBL, Qbl,4>bl) to harmonic coordinates {tn,XH, 
Uh,zh) was found by Cook and Scheel [p8| to be 



t'H = tBL 



■In 
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r + — r_ 

xh + iyn = {tbl —M + ia)e l(l> sm9 B L, 
zh = {tbl - M) cos 9 bl, 



where 
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(Bib) 
(Blc) 



(B2) 



M and a are the Kerr mass and spin parameters, and r± = M ± \J M 2 



The angular frequency in Boyer-Lindquist 



coordinates is the same as in harmonic coordinates, u> = d<j)BL I 'dt-BL ~ d(pfj/dtu, where (f)jj = tan 1 (jjh/xh)- The 
transformation of the radial coordinate can be written as 



r H = r B L - M + 0(a 2 ), 



(B3) 



which is sufficiently accurate for spin-orbit effects. 

The orbital angular frequency u> of a test particle in a circular orbit in the Kerr spacetime with spin perpendicular 
to the plane of the orbit is given by Eq. (4.26) of Ref. [p7| . That equation can be written as 
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Here q = a/M , sj_ is the magnitude of the spin of the particle divided by the square of its mass, m p is the particle's 
mass, and Mbh is the (much greater) black hole's mass. Note that the effect of the particle's spin s± is already 
0{m v IM.BH) smaller than that of the black hole.^ The parameter s± is explicitly written in terms of the coordinate 
components of the spin vector as 



S; 



BL + a COS#BL, 



where S® is the #BL-component of the spin vector of particle. Since we consider the equatorial plane, 
the coordinate transformation (Bl), we find that 



(B5) 
7r/2. Using 



4 The coordinate transformation (B3) is modified by the presence of the particle. However the modification appears at 
0{m p / Mbh) in the spinless terms and at 0(m p /MBH) 2 in the spin terms; thus it does not affect the comparison here. 
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(B6) 



where S% is the z-component of the spin vector of the particle in harmonic coordinate and m 2 = m p is the particle's 
mass. 

Using above formulas, we write the angular frequency of a test particle in harmonic coordinates as 
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where s 2 = S 2 /m\. By taking our post-Newtonian result ( |5.22j ) for u> 2 , setting body 1 to be the Kerr black hole and 
body 2 to be the particle, taking the test-mass limit 
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m 2 
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M, 



BH 

Xs = ^(q + s 2 ), 
Xa = ^{q~ s 2 ) , 

we obtain the same result. Thus our post-Newtonian calculation agrees (in the test-mass limit) with the previous 
results of black-hole perturbation theory. 



APPENDIX C: LORENTZ INVARIANCE OF THE EQUATIONS OF MOTION 



In this appendix, we briefly describe the Lorentz invariance of the equations of motion (5.5). It is well known 
that the 2PN equations of motion, including spin-orbit and spin-spin terms, are Lorentz invariant Therefore we 
concentrate on the new 2.5PN spin-orbit terms. 

Consider the case where the coordinates (t' , x % ) are moving relative to (t, x 1 ) with constant velocity U l . The Lorentz 
transform of a vector A is defined as 



A 11 = K» V ,A V 



(CI) 



where 



A - 
A\- 



o 



7, 

AV 

S ij - 



(7-1), 



(C2a) 
(C2b) 

(C2c) 



and 7 = 1/vl — U 2 . We expand the Lorentz transform (CI) in powers of U — \U % \, assuming that U is of the same 
order as va so that the slow-motion assumption underlying the post-Newtonian expansion is still valid. 

We also account for the fact that the time slices in the primed coordinates intersect the bodies' world lines at 
different points than in the unprimed coordinates. Without loss of generality, we can take the time slice in the 
unprimed coordinates to be t = and the primed time slice to be t' = 0. Then it follows from the inverse Lorentz 
transformation that quantities evaluated at body A at t = are evaluated at 



t' = t / A = y A (t = Q) = -U k y A (t = Q) 



(C3) 



When transforming the equations of motion to the primed coordinates, we must evaluate the primed vectors at t' = 
rather than t' — t' A . Therefore we expand the transformed equations of motion in powers of the small quantity t' A . 
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Specifically, the position y\ of body 1 on the time slice t — is related to the position y\ of body 1 on the slice 
t' = by 



y\^y\- -UWy\ - v\ Wy{ , 



(C4) 



where v\ = dy\ j dt' . Here and below, (unprimed) vectors on the left-hand sides of equations are evaluated at 



t = while (primed) vectors on right-hand sides are evaluated at t' 
four-acceleration, we find that 



= 0. By transforming the four- velocity and 
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KP u (—U)yi. In the same way, by considering the Lorentz 



transform and time slicing of the spin variables and the spin supplementary condition, we find that 



si 
si j 



st + ^ius'j + w'is'A), 



(C7) 



si- 



(v[U) 



jk v1 (S[U) 



]k u k [(v' 1 s' 1 ) + (us' 1 )}~u^'s( ]k 'u k 
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(C8) 



The transformations of the variables of body 2 are given by exchanging subscripts 1 and 2 in the above formulas. 

We insert these formulas into the equations of motion (5.5)— (5. 1C), keeping the lowest-order acceleration on the 
right hand side of Eq. (C6) but eliminating other (post-Newtonian) accelerations with the 1.5PN equations of motion 
in the primed coordina tes. T hus we obtain equations of motion in the primed coordinates which, to 2.5PN order, have 
the same form as Eq. (5.10). This completes the check of the Lorentz invariance of our 2.5PN equations of motion. 
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